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Abstract 

Assuming the existence of a fundamental thermodynamic relation, the classical thermodynamics 
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I. INTRODUCTION 


Classical macroscopic thermodynamics is a branch of physics of great generality, applica¬ 
ble to systems of elaborate structure, with different physical properties, such as mechanical, 
thermal, electric and magnetic. Therefore, thermodynamics possess great importance in 
many fields of physics and technical sciences 00 The task of thermodynamics is to define 
suitable macroscopic quantities to describe average properties of systems and explain how 
are these quantities related by means of universally valid equations. 

Even though thermodynamics is entirely of phenomenological nature, it is one of the 
best logically structured branches of physics. However, it is known, that systems which 
interact via the so-called long-range interactions, present non-conventional thermodynamic 
properties M. Gravitational, coulombic and dipolar systems, are some examples of long- 
range interactions, in which it is possible to observe non-common thermodynamic features, 
such as non-additivity, non-extensivity and the possibility of negative heat capacity. This 
work is inspired within this phenomenological thermodynamic spirit. 

Bekenstein pm proposed that an appropriate multiple of the event horizon area must 
be interpreted as the entropy. 




( 1 ) 


where S^h is the entropy of the black hole, ks is the Boltzmann constant, c is the speed of light 
in vacuum, h is the Planck’s constant divided by 2%, G is the universal constant of gravity 
and A is the area of the event horizon. Moreover, Hawking et al |B] proved that certain laws of 
black hole mechanics are mathematical analogous of zeroth and first laws of thermodynamics. 
Shortly after, it was proved that in addition to the laws of equilibrium thermodynamics, black 
holes obey as well, the standard theory of non-equilibrium thermodynamics in the form of a 
fluctuation-dissipation theorem P . Development of the so-called black hole thermodynamics 
still continues [THI - fll] . 

Existence of non-zero temperature for a black hole supports the hypothesis that the 
relation between entropy and the area of the event horizon of a black hole has some physical 
significance. Knowledge of this quantity makes possible to obtain the entropy of a black 
hole as a function of some fundamental parameters which describes such system; namely, a 
thermodynamic fundamental equation for the black hole. Thermodynamics is used to find 
information about a gravitational system composed by a black hole, approaching to the 
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problem from the classical thermodynamics point of view [I]. 

Although a formal study of black hole thermodynamics can be explored from a more 
fundamental perspective it is not the goal of this work to pursue this path. Instead, 

in this work exploration of the thermodynamic information that can be extracted from the 
assumed existence of a thermodynamic fundamental equation is considered, as would be done 
for an ordinary thermodynamic system. This approach proves to be successful to recover 
the mathematical structure, completely analogous, to conventional systems, in particular 
magnetic ones. 

This paper is organized as follows. In section two, the fundamental thermodynamic 
equation of a black hole, its natural parameters and the derived equations of state are briefly 
discussed. The non-common features observed in these equations are examined, including 
how can they be interpreted. Mathematical analogous of the Euler equation and the Gibbs- 
Duhem relation for these systems are calculated. In section three new response functions 
are defined and analyzed, in particular, specific heat and its odd behavior is studied. With 
the aid of TdS equations, independence of the new response functions are found. In section 
four, the three possible Legendre transformations are performed and the Maxwell relations 
are given. The mnemonic diagrams for those relations are constructed (thermodynamic 
squares) exhibiting an interesting relation with magnetic systems. A particular application 
of the formalism to the “rotocaloric” effect is discussed. In section five, stability conditions 
are reviewed and some consequences are analyzed. Some conclusions are given in section 
six. 

II. THE FUNDAMENTAL EQUATION: GIBBS-DUHEM AND EULER RELA¬ 
TIONS 

In this work, the simplest black holes, i.e., those with spherical symmetry |20] are con¬ 
sidered. The most general one is the Kerr-Newman black hole, which is characterized by 
three parameters: mass M, angular momentum J, and electric charge Q. Fundamental 
equation in energy representation can be found by substituting the area A(M, J, Q) in Q, 
using U = Mc^, and inverting the expression to get U = U{S, J, Q): 
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The first-order differential for the internal energy is: 


dU 





( 3 ) 


In analogy with internal energy of a thermodynamic simple system, U = U{S,V, Ni), 
where V stands for volume and Ni for mole number of the z-th component, partial derivatives 
appearing in ([^ play the role of intensive parameters |0|, 


T = 





(4) 


Where T is the temperature of the black hole, its angular velocity and $ describes the 
electric potential for a charge Q evaluated at a distance equal to the event horizon radius [21] . 
These relations are equivalent to the EOS in a simple thermodynamic system. Then, we 
have the expression: 


dU = TdS + ndJ + ^dQ. 


(5) 


Such relation is analogous to the thermodynamic fundamental equation of a simple sys¬ 
tem. In a conventional thermodynamic system, derivatives of internal energy with respect 
to extensive parameters give intensive parameters. This property is reflected in the fact 
that for a conventional thermodynamic system, the fundamental equation is a homogeneous 
first-order function of some extensive parameters. Similarly, the EOS of a conventional ther¬ 
modynamic system are homogeneous zero order functions. However, for the Kerr-Newman 
black hole, fundamental equation Q, is not an homogeneous first-order function. Two major 
implications for thermodynamics systems: first, it is an indication of non-additivity; second, 
non-extensivity is present, that implies that the EOS are not homogeneous zero-order func¬ 
tions. 

The fundamental cause of those particular thermodynamic properties lies in the gravita¬ 
tional nature of the system. Gravitation falls in the so-called long-range interactions. It is 
largely known that thermodynamical description of such systems presents several complica¬ 
tions, among them, non-additivity and broken extensivity [3HS]. 

Eor the Kerr-Newman black hole, EOS obtained by derivation of eq. are: 


n 


f J \ 

y hG s) 


• E(5,j,g), 


( 6 ) 
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( 7 ) 


/ T^ksC^Q TT‘^k%C^Q^\ 

y G ^ hGS J 


T = 



E(S,J,Q) 


(?{hcS + TrksQ"^) 


with S(5', J, Q) given by: 


c 


[{hcS + nksQ^Y + (SttAibcJ)^] 


1 


S = 



(hcS + TvksQ'^Y + (27rfcscJ)^ 



( 8 ) 


(9) 


Given the astrophysical relevance of the angular momentum compared to the electric charge, 
henceforth in this work a Kerr black hole is considered, characterized by its mass M, and 
angular momentum J. As a consequence, a fundamental equation U = U{S, J) by taking 
Q = 0 in eq. (|^ is used, for the sake of simplicity. 

We analyze (in natural units) the EOS, T = T{S,J) and fl = fl{S,J) by plotting 
each one with a fixed variable for each curve and studying its behavior. Starting with the 
angular velocity fl(S, J), in FIG. 1, it is shown that for fixed entropy and increasing angular 
momentum, the angular velocity increases asymptotically until it reaches a maximum value. 
The rise of this asymptote can be understood by the physical limit imposed to a rotating 
black hole, which will be discussed later. 

Gonsidering FIG. 2, if the angular momentum is kept constant for an increasing entropy, 
the angular velocity decreases to zero asymptotically for maximum entropy. The system 
can produce large amounts of energy at the expense of its rotational energy, maximizing 
entropy [22]. 

Similarly, for the temperature T{S, J), there are two cases: In FIG. 3, if S is kept constant it 
is possible to notice that, for an increasing angular momentum, temperature falls rapidly to 
zero, therefore, the angular momentum plays the role of a temperature attenuator. During 
the process of energy production at expense of rotational energy, this energy is radiated as 
thermic radiation, reducing J and increasing T. 

Considering FIG. 4, for a fixed J, temperature reaches a maximum as entropy increases, 
in this plot there are two regions due to the competition between J and S. There is a 
region where temperature is a monotonically increasing function of entropy, and a second 
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region where T drops asymptotically until it reaches a certain value for maximum entropy. 
Presence of those regions for temperature is related to the stability of the system. 

In the entropic representation, the system has a relation between internal energy and 
angular momentum, establishing a boundary beyond which the proposed model has no 
physical meaning. That happens since reality condition for entropy is violated in this limit, 

Gc-^U^ > J. (10) 


This expression can be interpreted saying that for a fixed value of the internal energy, angular 


momentum has a limiting value, given by eq. (10). Behavior of the angular momentum and 
the allowed region for internal energy are shown in FIG. 5. It is important to mention 
that this maximum value of J also appears in general relativity, referring to extreme Kerr 
black holes and naked singularities [23]. Existence of the same prediction in two different 
approaches is a good indication that postulation of a fundamental thermodynamic equation 
reasonably captures the essence of the system. 

Following the canonical thermodynamic formalism, analogous of Euler equation, since ex- 
tensivity is broken, is calculated. The internal energy is an homogeneous 1/2-order function, 
which implies that: 

2TS + 2VLJ ^U. (11) 

From this expression, a Gibbs-Duhem-like relation can be calculated. 


dU ^ - 2 {SdT + Jdny, 


( 12 ) 


with this procedure, a relation for energy as a function of temperature and angular veloc¬ 
ity 17 (T, Q) is obtained. However, eq. 0 it is not a fundamental thermodynamic 

equation, for this function it follows: 


dU 

&T 


-2S] 



(13) 


III. RESPONSE FUNCTIONS AND TdS EQUATIONS 

Second derivatives of the fundamental equation, also called response functions, are usu¬ 
ally descriptive of material properties and often these quantities are of the most direct 
physical interest. For the Kerr black hole we proposed the following basic set of response 
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functions, composed by olq^xt^Cq. With this particular set, it is possible to recover all 
the thermodynamic information. As we will show, dehnitions given below reproduces the 
conventional structure of thermodynamics, showing a close analogy with the thermodynam¬ 
ics of magnetic systems but different to the thermodynamics of fluids. To our knowledge, 
description of thermodynamics for the Kerr black hole, using this set of response functions, 
it is presented for the first time. 

Coefficient is defined by: 

The quantity is the change in the angular momentum per unit change in the temperature 
of a Kerr black hole maintained at constant angular velocity. 

The response function xt is given by: 


Xt — 


dJ 


(15) 


\ / T 

Xt is the increase in the angular momentum per unit change in the angular velocity at 
constant temperature, and can be interpreted as an isothermic rotational susceptibility. 
The heat capacity at constant angular velocity is defined by: 


Cn = T 


dS 

WT 


dQ 

dT 


(16) 


\ /a \ /a 

The heat capacity at constant angular velocity is the quasi-static heat flux required to 
produce unit change in the temperature of a Kerr black hole maintained at constant angular 
velocity. 

All the other response functions of the Kerr black hole can be expressed in terms of this 
basic set Cq,Xt,cx.q.- This set of functions are in fact, the three possible second derivatives 
of the thermodynamic potential G = G{T^ Vt) that will be introduced in section three, 


Xt — — 



( d^G\ 

“ [WM j ’ 



(17) 


Two additional response functions of interest are the heat capacity at constant angular 


momentum Cj, and the isentropic susceptibility, xs- 


The heat capacity at constant angular momentum, defined by: 
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J 


(18) 


C, = T 




I dT 




dT 


it is a measure of the quasi-static heat flux needed to change in one unit the temperature of 
a Kerr black hole maintained at constant angular momentum. 

The response function xs is given by: 


(19) 

s 

This quantity characterizes the change in the angular momentum associated with an isen- 
tropic change in the angular velocity. 

The TdS equations are necessary to obtain relations between response functions. These 
equations can be found by considering entropy as a function of certain natural variables, 
differentiating and rearranging the expression to find a relation between TdS and certain 
second derivatives. For a Kerr black hole the first TdS equation can be found by considering 
entropy as a function of the temperature and the angular velocity, S = S(T,fl), 



TdS 


CndT + T 





( 20 ) 


\ /a 

the second TdS equation can be obtained considering entropy as a function of the temper¬ 
ature and the angular momentum, S = S{T,J), 


TdS = CjdT 


T 


on 


dJ. 


( 21 ) 


\ / J 

Relations between response functions can be found through manipulation of TdS equations 
and the use of Maxwell relations, discussed in section four. Two relations that describe the 
dependence of Cj and xs in terms of the basic set of response functions are: 


Xt{Cq, — Cj) — TaQ, 

CnixT — Xs) = ToiQ- 

The specific heat at constant angular momentum is calculated, 


Cr = 


du 

WT 


hc^ 


2T^ \ ATTkeG HGS^ 


( 22 ) 

(23) 


J 


(24) 



It is possible to notice the existence of a negative region. Since the second term is reduced 
for an increasing entropy and subtraction inside the brackets will eventually result in a posi¬ 
tive quantity, causing that the specific heat becomes negative, for sufficiently large values of 
entropy. Negative region for the specific heat is known long ago for gravitational systems |24] 
in general, and for black holes in particular ISlElIiniES]. Two interpretations about negative 
specific heat arise. 

One is related to the transfer of energy as heat; having a region of negative heat implies 
that the temperature of a black hole raises to a higher value, with energy transfer [2H]- 
The second one is related to thermodynamic stability of the system. For a negative spe¬ 
cific heat, thermodynamic systems are in unstable states passing perhaps, through a phase 
transition |28j . 

IV. THERMODYNAMIC POTENTIALS AND MAXWELL RELATIONS 

Another point of interest are the thermodynamic potentials. Alternative representations 
for the thermodynamic fundamental equation, can be obtained via the application of a 
Legendre transformation to U = U{S,J). For a Kerr black hole there are three possible 
Legendre transformation. 

Starting with the mathematical analogous of the Helmholtz free energy, F = F{T, J) = 
U[T], 


F = U -TS, 


and 


dT' 


(25) 


which can be written as: 


F = S(A, J) 


hF 


l^ + 


3J2 


hc^ 


( 26 ) 


AirksG I 2S \ 27rA:B 
with S(S', J) given by eq. i§. for (5 = 0. It is worth to mention that eq. (|2^ must be a 


function of temperature and angular momentum as independent variables; instead, eq. (26) 


is a function of the entropy and the angular momentum, F = F{S, J). We were not able to 
find an analytical expression for Helmholtz free energy in that terms. 

The Helmholtz free energy for a Kerr black hole in diathermal contact with a heat reservoir 
can be interpreted as the available work at a constant temperature, since dW = fldJ. 
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A second Legendre transformation, can be obtained replacing the angular momentum 
J, by the angular velocity fl, as the independent variable in fundamental equation, H = 
H{S, fi) = U[fl], this leads to: 


H = u-nj, 

the potential H is given by: 


and 


- J = 


H = 




1 + 


^2 


hGS 


- 0 ? 


OH 


1/2 


(27) 


(28) 


which is a function of entropy and angular velocity as independent variables. Therefore, eq. 


(28) is an explicit expression for the potential H = H{S,n). It can be interpreted as the 
heat added to a Kerr black hole at constant angular velocity, since dQ = TdS. 

The last thermodynamic potential is obtained by simultaneously replacing both, the 
entropy and the angular momentum, by the temperature and the angular velocity as inde¬ 
pendent variables in the fundamental equation, G = G{T, fl) = U[T, fl]. 


G = U -TS -nj, and 
an expression for the G potential is given by: 


_ T-— -Q-^. 

“ “ M’ ~ ~&T'' 


(29) 


G = 


1 


TrkBC^ \ ( ^2 , t 2 2 

HGS J \^47r2fc| ^ 


1 1/2 


(30) 


As for Helmholtz free energy in eq. (26), an expression for G(T, Q) as a function of S and 
J instead of its natural variables was obtained. Physical interpretation of this potential 
for a Kerr black hole it is not clear, since in a conventional thermodynamic system, the 
mathematical analogous of this potential used to be associated with a third extensive pa¬ 
rameter, missing in this case. However, for a Kerr-Newman black hole, this third parameter 
is the electric charge Q, therefore, this potential can be interpreted as the available work 
at constant temperature and angular velocity, since dW = d>d(5, for a Kerr-Newman black 
hole. 

Different representations for the fundamental equation leads to the existence of various 
derivatives of independent variables related to each thermodynamic potential, apparently 
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unrelated with other derivatives. However, there are some relations among such deriva¬ 
tives; hence, just a few of them can be considered as independent. These relations between 
derivatives are known as Maxwell relations and arise from the equality of the mixed partial 
derivatives of the fundamental equation, expressed in any of its alternative representations. 
Since there are only two independent variables, no more than a pair of mixed derivatives 
exists in any representation of the fundamental equation. For the energetic representation 
U = U{S,J), Maxwell relations are: 


{s,n) 

{T,n) 

{T,J) 

{S,J) 


d^H 
dSdn ~~ 

d^G 

dTdil 

d^F 
dJdT '' 

d‘^U 



(31) 

(32) 

(33) 

(34) 


dJdS dSdJ 

Maxwell relations can be remembered conveniently in terms of simple mnemonic dia¬ 
grams [29]. In this case, the construction rules of those diagrams are inspired in the con¬ 
struction of the same mnemonic diagrams for magnetic systems HU: A minus sign is added 
to the extensive variable of magnetic work, with this adaptation it is possible to build the 
mnemonic diagrams for magnetic systems. For Kerr black holes, the same result can be 
achieved by adding a minus sign to the angular momentum. Maxwell relations are repre¬ 


sented in the next diagrams. In FIG. 7, the Maxwell relation for the potential H, (31) is 
depicted. The remaining Maxwell relations are represented in FIG. 8 - FIG. 10 as successive 
rotations of the thermodynamic square shown in FIG. 7. 


A. The “Rotocaloric Effect” 

An additional example of the similarities between thermodynamics of magnetic systems 
and thermodynamics of black holes, lies in the exploration of the magnetocaloric effect 
(MCE) and its counterpart for a Kerr black hole, here named as “rotocaloric effect”. 

The MCE is a magnetothermal phenomena exhibited by magnetic materials, in which a 
reversible change in temperature of the material occurs under adiabatic and usually isobaric 
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conditions in the presence of a variable magnetic field EH For a magnetic system with a 
total differential for entropy expressed as a function of temperature T, and the magnetic 
field H, the change in temperature under variations of H, can be obtained from: 


dT _ {dS/dH)T 
m ~ ~ {dS/dT)H ■ 


For a Kerr black hole a similar process could appear, a reversible change in temperature 
of the event horizon under adiabatic conditions due to variations in the angular momentum 


of the black hole. To calculate this changes in T, the TdS equation (21), is used. In an 
adiabatic process TdS = 0: 




(36) 


using the definition of Cj and the Maxwell relations, it is possible to express the variation 
in T as: 


dT ^ {dS/dJ)T . . 

dJ {dS/dT)j' ^ ’ 

To compare both phenomena, the MCE is applied to a paramagnetic substance composed 

by N freely removable dipoles with a permanent magnetic dipole moment fl [2]. Consid¬ 
ering only the dominant term in entropy S{T, H), the variation in the temperature of this 
paramagnetic substance under changes in the applied magnetic field, is given by: 


dT _T 
dH~H' 


(38) 


Similarly, to calculate the value of the rotocaloric effect for a Kerr black hole, an explicit 
relation for entropy as a function of temperature and angular momentum is required. The 
equation of state for temperature eq. (|^, considering <5 = 0, was used to find this relation. 
It is necessary to invert T = T{S,J) to obtain S = S{T,J). An approximated relation 
S{T, J), was estimated considering an iterative process. 

For a Kerr black hole, the first-order approximation was calculated proposing Sq as the 
entropy of a Schwarzschild black hole, the simplest one, by considering J = Q = 0 in eq. 
([^, and taking the derivative with respect of temperature: 


SoiT) 


hc^ 1 
IQi^GhsT^' 


(39) 


12 






With this relation, the first-order approximation <S'i(r, J), for the exact relation S{T, J), is 
given by: 


.5i = 


7i 


+ 


72 


( 40 ) 


J47'10 ' J2j’e’ 

where 71 and 72 are a combination of the fundamental constants used in this work. If only 


the dominant term 1/J^T^ is considered in eq. (37), the variation in the temperature under 
changes in the angular momentum for a Kerr black hole is: 


(41) 


dT _ T 

Jj ~ 

As J and T are non-negative in all the entropy domain, the rotocaloric effect produces a 
drop in the temperature of the black hole for an increase in the angular momentum and vice 
versa, this result was previously obtained analyzing the equation of state for the angular 
velocity, eq. ([^ in section two. The results in eqs. (38) and (41) are very similar, but 
with different sign. Therefore the changes in temperature work in the opposite direction for 
rotocaloric effect, compared with magnetic systems. 

Due the analogies between MCE and rotocaloric effect, it is possible to construct a re¬ 
frigeration cycle for the Kerr black hole. This process is very similar than the magnetic one, 
only with the opposite sign. 

A deeper exploration of these effects will be given in a subsequent work. 


V. STABILITY CONDITIONS 

As is discussed in section three, exist a negative region for the specific heat Cj, which 
indicates non-stability in the system. It is convenient to use the entropic representation of 
the thermodynamic fundamental equation to analyze stability conditions. 


A = 


27rfci 

he 


Gc-^U^ + ( Gc-^U"^ - c^J 


-87-r4 t2 


1/2 


The system is thermodynamically stable one, if the entropy S{U, J) satisfies: 


( 42 ) 


S(U, J + AJ) + S(U, J -AJ)< 2S(U, J), 
S(U + AU, J) + S{U - AU, J) < 2S(U, J); 


( 43 ) 

( 44 ) 
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or for infinitesimal displacements: 


d^S 


< 0 , 


d^S 


< 0. 


(46) 


\ / J \ / u 

For a stable system, the graphical representation of the previous conditions results in a 
concave curve for entropy. If the entropy has a non-concave curve, the related system is 
unstable. When the entropy possess both, concave and non-concave regions, the system is 
not stable; instead, it is a locally stable one, but globally unstable. 

Stability under variations in the angular momentum, at a fixed internal energy, is analyzed 
in FIG. 11. Entropy is a concave function in all its range of definition for the exchange of 
angular momentum at constant internal energy. Therefore, the Kerr black hole is stable 
under changes in the angular momentum. 

If the angular momentum is kept constant and variations of the internal energy are allowed, 
it is possible to notice that the entropy is non-concave in almost all its range, except for a 
region of very low energy, where the entropy is a concave function. Hence, a Kerr black hole 
is locally stable around K ~ 0, but it is globally unstable under variations in its internal 
energy, as can be seen in FIG. 12. 

If variations in both, the angular momentum and the internal energy are allowed, the 
three-dimensional space (5', G, J), must be considered. The global stability condition requires 
that the entropy surface S{U, J) lies everywhere below its tangent planes. Accordingly, for 
arbitrary AG and AJ: 


Siu + AG, J + AJ) + S{U -AU,J- AJ) < 2S{U, J); (46) 

and locally, for AG —)-0AG ^0: 

d^S d^S f ^ ^ . 

~ ydUdJj - 

Behavior of the entropy in the subspace {S, G, J) for a Kerr black hole, is shown in FIG. 13. 
Entropy is a non-concave function of G and J; in correspondence with two-dimensional case 
for the entropy as a function of the internal energy, a region of stability exists. This locally 
stable region is available for small values of energy and angular momentum. 

Evidence that a Kerr black hole is a thermodynamic system in unstable equilibrium in 
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general, but in a state of locally stable equilibrium for very small values of its mass and 
angular momentum, is founded. 

VI. CONCLUSIONS 

The implementation of the conventional thermodynamic formalism for black holes is 
permitted by the postulation of a fundamental thermodynamic equation of the form 
U = U{S,J,Q). The proposal of considering a purely thermodynamic vision appears to 
be fruitful. The application of this formalism in such systems presents particular complica¬ 
tions, as non-additivity or instability of equilibrium states, typical of a system with dominant 
gravitational interactions; however, this development shows the presence of considerations 
suggesting that, beyond the mathematical exercise, the proposed model really represents 
some aspects of a very interesting behavior. This can be seen in the predictions made, that 
can be found in other frameworks of physics and in resemble with the thermodynamics of 
magnetic systems. The presented formalism could be a start point for the research of inter¬ 
esting analogous thermodynamic processes for black holes, as the Joule-Thomson effect and 
others. Results about these topics will be presented in a future work. 
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Figure captions 


Figure 1. Angular velocity vs. angular momentum for a given entropy. Asymptote in 
fl is linked to the so-called Extreme Kerr black hole. 

Figure 2. Angular velocity at a given J. Rotation rate decreases as the thermal energy 
is radiated increasing entropy. 

Figure 3. Temperature at fixed S. Inverse function J = J(T) is plotted to show that 
T increases when the rotating black hole transforms its rotational energy into thermal 
radiation. 


Figure 4. Temperature at given J. The maximum is related with stability and the sign 
of the specific heat. 


Figure 5. Forbidden region (over the curve) for angular momentum. 


Figure 6. Specific heat Cj showing negative values for large enough entropy. 


Figure 7. Thermodynamic square representing Maxwell relation (31) for the thermo¬ 
dynamic potential H = H(S,fl). 


Figure 8 . Thermodynamic square representing Maxwell relation (32) for the thermo¬ 
dynamic potential G = G{T, fl). 


Figure 9. Thermodynamic square of Maxwell relation (33) for Helmholtz free energy, 
F = F(T, J). 


Figure 10. Thermodynamic square representation of Maxwell relation (34) for internal 
energy, U = U{S,J). 


Figure 11. Entropic representation of fundamental equation S{U, J) at constant energy. 
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Concavity of the curve implies stability. 


Figure 12. Entropy S{U,J) at a fixed J. Entropy is non-concave at almost all their 
domain, except for a small region of low energy. 

Figure 13. Three-dimensional space {S,U,J), showing unstable equilibrium for the 
system. 
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